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FOREWORD 


This  report  was  prepared  under  Air  Force  Contract  No.  AF  08(635)-171S, 
(U)  "Theory  of  Speed  Impact).  The  work  was  administered  under  the 
direction  of  APGC  (PGTWR).  The  project  engineer  was  Dr.  T.  D.  Riney. 


The  assistance  of  Mr.  P.  R.  Chemoff  of  the  General  Electric  Co.  in  the 
preparation  of  this  report  is  gratefully  acknowledged.  He  contributed  to  the 
numerical  analysis  and,  also,  programed  the  dUCsrenoe  equations  for  machine 
calculations. 


ABSTRACT 


I 


An  explicit  finite  difference  formulation  of  the  equations  governing  the 
one -dimensional  visco-plastic  model  was  presented  in  the  Second  Quarterly 
Report.  The  scheme  has  been  programed  on  the  IBM  7090  and  e3q)loratory  cal¬ 
culations  made  for  several  values  of  the  parameters  fio  (viscosity  factor),  r  o 
(yield  stress)  and  vq  ^nqmct  velocity).  Excellent  refindts  are  obtained  for  cer¬ 
tain  ranges  of  the  parameters:  these  results  are  discussed  in  detail.  Tw  larger 
valties  of  fio  and  v^,  however,  tiie  restriction  on  the  size  of  At  is  very  severe. 

A  stidbiUty  axialysls  is  performed  which  explains  this  requirement. 

To  circumvent  this  difficulty  the  governing  equations  are  recast  in  an 
implicit  finite  difference  scheme  which  is  only  ^id  for  Mo  >  viiich  is  then 

unconditionally  stable.  The  method  is  described  in  detail  and  a  flow  chart  is 
given  fbr  tbe  program  vddoh  is  currently  being  written. 
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K  number  of  iterations  performed 

J  position  station  chosen  to  be  in  front  of  disturbance 
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8  quantity  defined  in  (A-16) 

5  quantity  defined  in  (A-16) 

G  amplification  matrix 

X  eigenvalues  of  amplification  matrix 


V 


TABLE  OF  CONTENTS 


Page 


INTRODUCTION .  1 

GENERAL  REMARKS. .  1 

(a)  Choice  of  Parameters .  1 

(b)  Description  of  the  Flow .  2 


EXPLICIT  DIFFERENCE  SCHEME . 4 

(a)  Stability  Criteria . 4 

(b)  Results  of  Calculations . y 

IMPLICIT  DIFFERENCE  SCHEME .  9 

(a)  Difference  Equations . o 


(b)  Initial  and  Symmetry  Conditions. . 

(c)  Iteration  Procedure . 

(d)  Two-Term  Recurrence  Relations. 

(e)  Three -Term  Recurrence  Relation 

CONCLUSIONS . . 


BIBLIOGRAPHY .  I9 

APPENDIX . . 


LIST  OF  ILLUSTRATIONS 


♦ 


i 


Figure 

1  Calculated  pressure  pro&le  (a=l.  5,  ^=0.1,  M=0.  05)  compared 

with  the  Rankine-Hugoniot  solution 

2  Calculated  and  Rankine-Hugoniot  values  of  the  shock  velocity 

relative  to  the  interface 

3  Pvilse  velocities  at  various  times  compared  with  asymptotic 

values 

4  Pulse  velocities  at  various  times  compared  with  asymptotic 

values 

5  Typical  plot  of  von  Mises  flow  statistic  at  various  positions 

of  the  disturbance 

6  Illustration  of  impact  situation  showing  (a)  center  of  mass 

coordinates,  (b)  Eulerian  space  -  mesh  points 

7  Tabular  display  of  initial  and  boundary  data  for  implicit 

difference  scheme 

8  Flow  chart  depicting  the  sequence  of  calculations  in  the 

implicit  difference  scheme 


vii 


INTRODUCTION 


In  the  First  Quarterly  Report  (Ref.  1)  a  mathematical  model  was  proposed 
for  the  study  of  hypervelocity  intact  which  takes  into  account  the  effects  of 
material  viscosity  and  dynamic  yield  strength.  In  this  formulation  the  noaterial 
is  considered  rigid  if  stressed  below  Its  flow  point,  whereas  the  material  acts 
like  a  Newtonian  viscous  liquid  when  stressed  above  this  value. 

In  the  Second  Quarterly  Report  (Ref.  2)  dimensionless  parameters  were 
found  which  control  the  relative  Importance  of  the  inertial,  viscous  and  strength 
terms  in  the  governing  partial  differential  equations.  Both  the  viscosity  factor 
Mo  aiid  the  dynamic  yield  stress  Tq  appear  in  these  parameters.  Since  definitive 
data  is  not  available  for  either  of  tiiese  quantities  a  proposal  was  made  that 
exploratory  calculations  be  performed  using  a  one-dimensional  model  in  which 
the  values  of  Mo  ^o  varied  over  several  orders  of  magnitude.  For  this 
purpose  the  governing  equations  were  written  in  one-dimensional  Lagrangian 
form  and  then  recast  in  finite  difference  forni  in  preparation  for  the  exploratory 
calculations.  In  this  scheme  the  dependent  variables  at  time  (n  i)  are  ex¬ 
pressed  explicitly  in  terms  of  their  values  at  time  nAt. 

The  present  report  opens  with  a  description  of  some  of  the  results  of  the 
calculations  made  on  the  IBM  7090.  Very  good  results  are  obtained  for  certain 
ranges  of  the  parameters  Mo«  tq  and  intact  velocity  Vq.  For  larger  values  of 
Mo  and  vq,  however,  the  calculations  do  not  converge  but  oscillate  with  Increas¬ 
ing  amplitude  as  n  Increases.  A  stability  analysis  is  given  which  shows  that 
this  effect  results  from  the  inclusion  of  the  viscosity  which  causes  the  usual 
restriction  on  the  ratio  At/Ax  to  be  replaced  by  a  restriction  on  the  ratio  on 
At/(Ax)^.  This  conclusion  is  verified  by  the  preliminary  calculations.  To  cir¬ 
cumvent  this  difficulty  the  governing  equations  are  written  in  one-dimensional 
Eulerian  form  and  then  recast  in  a  finite  scheme  which  is  only  valid  when  Hfy>0, 
but  which  is  then  unconditionally  stable.  Here  the  dependent  variables  enter 
inqjlicitly  so  that  an  iterative  procedure  is  necessary  in  proceeding  from  time 
n  to  time  n  +  1.  The  method  is  described  in  detail  and  a  flow  chart  is  given  for 
the  calculations. 


GENERAL  REMARKS 


(a)  Choice  of  Parameters 


In  choosing  trial  values  for  Mq  starting  point  has  been  the  values  quoted 
by  Perzyna  (Ref.  3) .  On  the  basis  of  the  results  of  internal  energy  measure¬ 
ments  by  Kolsky  (Ref.  4)  he  computes  the  approximate  values^ 


^The  gram-centimeter- microsecond  system  is  used  throughout  this  report. 
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steel:  =  0. 8  gm  cm“^  microsec" 

Copper:  fiQ  =  0.4  gm  cm“^  microsec”^. 

CXir  calculations  have  all  been  for  iron  and  Hq  is  assumed  to  be  within  a  factor 
of  ten  of  0 . 8 . 

For  mild  steel  Reiner  (Ref.  5)  quotes  the  static  yield  stress  to  be  approxi¬ 
mately 

o  1  —2 

Tq  =  10  gm  cm  ^  mlcrosec.  , 

i.e.  ten  kilobars.  For  our  calculations  it  is  the  dynamic  yield  stress  that  is 
pertinent.  The  above  value  is  assumed  to  be  the  lowest  value  likely,  aikd  Tq  is 
varied  up  to  one  hundred  times  as  large,  i.e.  one  megabar. 


The  various  combinations  of  assumptions  for  Tq  and  impact  velocity, 

Vq,  are  displayed  in  Table  I.  The  choices  of  Vq  (  =  0.  5,  4,  7.  5  cm/microsec) 
represent  thf  extremes  and  mean  of  the  meteoroid  velocity  range.  The  values 
of  the  dimen.  mless  paranubers  Bq  and  Rq  are  also  listed  for  each  combination 
of  1^,  Tq,  Vq.  They  are  seen  to  vary  widely  with  the  choice  of  combinations. 

(b)  Description  of  the  Flow 

If  the  two  impacting  bodies  are  perfect  fluids,  shocks  are  produced  which 
propagate  at  a  constant  velocity  and  the  jump  discontinuities  across  the  shocks 
remain  fixed.  With  viscosity  present  energy  is  being  dissipated  and  it  might 
be  expected  that  the  disturbance  would  be  attenuated  as  it  propagates .  In  the 
case  of  one-dimensional  impact  between  semi-infinite  bodies,  however,  the 
dissipation  effect  is  counteracted.  Here  there  is  an  infinite  reservoir  of  energy, 
a  steady  state  obtains  in  which  the  rate  of  increase  of  internal  energy  behind 
the  disturbance  just  equals  the  kinetic  energy  of  the  medium  entering  the  dis¬ 
turbance.  Consequently,  the  shape,  magnitude  and  velocity  of  the  disturbance 
(or  "pressure  pulse")  does  not  change  as  it  propagates;  These  characteristics 
of  the  steady-state  disturbance,  however,  are  affected  by  the  values  of  the 
parameters  Pq,  Tq  and  Vq. 

The  viscosity  coefficient  and  von  Mises  flow  rule  are  (Ref.  2,  p.  12) 


H- 


P 


dq| 

dxl 


*  Some  time,  of  course,  is  required 
established. 


for  a  stable  profile  to  be 
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TABLE  I 


Enumeration  of  parameter  coml 
microsecond  system,  p  =  7.  8 
a  is  chosen  as  1.  S.  -C.  =  1.  5  Ax. 


parameter  combinations  considered  for  iron-iron  impact.  The  units  are  in  the  gram-centimeter 
stem,  p  =  7.  8  gm/cm'^,  D  si.  In  computing  the  required  net  size  for  the  explicit  schemei 


Bingham  No 


Reynolds  No. 


Stability  Criteria 


Net  Size 


and  hence 


2  >  2 


,2.4.  Z  P 


(ii)- 


ro  +;f  H  ^0 


M 

dx 


*  T  (x)  (If-)' 


In  front  of  the  disturbance  3q/3x=0,  in  the  disturbance  5q/(3x<  0  and,  finally, 
3q/dx  drops  back  to  zero  after  the  disturbance  oasses.  Consequently, 
t2>  Tq^  in  the  disturbance,  but  drops  to  To^  after  it  passes.  This  means, 
in  terms  of  our  visco-plastic  model,  that  there  is  flow  only  in  that  part  of  the 
medium  through  which  the  disturbance  is  currently  passing;  it  again  becomes 
rigid  behind  the  disturbance. 

EXPLICIT  DIFFERENCE  SCHEME 


(a)  Stability  Criteria. 


In  the  Second  Quarterly  Report  a  trial  stability  criterion  was  suggested 
for  the  finite  difference  scheme  presented  there.  It  was  deduced  by  analogy 
from  a  criterion  derived  for  the  perfect  fluid  equations.  In  performing  pre¬ 
liminary  machine  calculations  the  trial  criterion  was  found  to  fail  for  the  higher 
values  of  the  impact  velocity  and  viscosity  parameters.  Therefore,  this  phe¬ 
nomenon  had  to  be  investigated  more  thoroughly.  In  the  Appendix  the  stability 
analysis  is  presented  in  some  detail;  the  essential  results  are  contained  in 
(A-22)  and  (A-24). 


To  apply  criterion  (A-22)  we  recall  that 

av  >  aq 

a  t  ’  ^0  ^  X 

and  approximate  b^/b'x.  by-VQ/2Ax.  Here  Vq/2  represents  the  velocity  of  the 
interface  as  calculated  from  the  Rankine-Hugoniot  conditions  (Ref.  2  eq.  64), 
In  our  calculations  we  have  chosenV^a  Ax,  a  =  1. 5,  so  that  the  criterion  re¬ 
duces  to 


I  /*o+4.5  Po  %  Ax 


(1) 


AL 

(Ax)^ 
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The  value  for  V  =  \/p  should  be  taken  from  the  Hugoniot  curve  for  the  particular 
material  (Ref.  2,  Fig.  4). 

Now,  if  we  are  away  from  the  shock  p  »  Q  +  n/^  and  condition  (A -24) 
becomes  roughly 

At  ^  ^0 

2  x/z  -  P  I 

^ ^ i^v  ^  aiii 

2 

The  adiabatic  sound  speed  c  is  equal  to 

Consequently,  in  regions  other  than  those  containing  a  shock, 


is  the  approximate  stability  criterion.  This  is  the  tentative  requirement  for 
stability  suggested  in  the  Second  Quarterly  Report  (Ref.  2,  eq.  66). 

If  (Vq/2)>c,  a  shock  will  occur  in  the  impacting  materials  and  criterion  (1) 
must  be  applied  even  when  Mq  =  0*  I**  the  case  of  iron,  for  example,  c  ^  0. 5  cm/ 
microsec.  and,  consequently,  criterion  (2)  holds  when  Vq  =  0.5,  while  at  impact 
velocities  4  and  7.5  criterion  (1)  obtains.  The  stability  criteria  have  been  cal¬ 
culated  for  each  of  the  parameter  combinations  considered.  The  results  are  listed 
in  Table  1  along  with  the  corresponding  value  of  At  for  the  choice  Ax  =  0. 1  cm. 


The  above  conclusions  have  been  substantiated  by  a  number  of  machine  cal¬ 
culations,  the  results  are  displayed  in  Table  II.  It  is  seen  that  the  calculated 
requirement  on  At  is  a  conservative  estimate  in  each  case.  It  may  be  expected 
that  the  theoretical  stability  criterion  is  close  to  the  required  condition  for  all  the 


2 

In  an  adiabatic  process  dU  +  p  dV  =  0.  Since  p  =  g  (V,  U  )  we  may  write 
dp  =  dVdg/dV  +  dU  dg/dV-  Elimination  of  dU  between  the  two  yields 
dp  =  (dg/dv  -  pdg/dv)  dv  =  v2  (pdg/du  -  dg/dv)  dp,  whence 

c^  =  dp/dp  =  y^ipdg/du-dg/d\). 
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TABLE  II 


Summary  of  preliminary  runs  testing  stability  of  explicit  difference  scheme 
applied  to  iron-iron  impact,  a  =  1.  5,  -t.  =  1.  5  Ax,  Ax  =  0.  1. 


NO. 

NET  SIZE 

RESULT 

1  - - - - 

Theory 

Trial 

1 

0.  06 

0.  100 

Unstable 

II- 

It 

0.  050 

Stable 

II 

It 

0.  065 

Stable 

2 

0.  027 

0.  05 

Stable 

3 

0.014 

0.  02 

Stable 

4 

0.0023 

0.  005 

Unstable 

It 

II 

0.  0025 

Unstable 

11 

II 

0.  002 

Stable 

13 

0.  06 

0.  05 

Stable 

17 

0. 0021 

0.  05 

Unstable 

If 

II 

0.  025 

Un  stable 

1 1 

1  1 

0.  01 

Unstable 

11 

M 

0.005 

Stable 

-4 

parameter  combinations.  For  the  most  severe  case,  At  =  5. 1  x  10  ,  this 

means  that  approximately  10,000  cycles  are  required  for  a  5  microsec.  run.  If 
70  space  mesh  points  are  required  this  means  700,000  point  calculations  which 
would  require  approximately  1.  4  hours  on  the  IBM  7090.  An  inordinate  amount  of 
machine  time  might  therefore  be  needed  to  complete  the  desired  production  runs. 
An  alternate  implicit  difference  scheme  is  therefore  proposed  for  those  parameter 
combinations  with  larger  Hq  and  Vq  valves.  Before  turning  to  this,  however,  let 
us  first  examine  the  results  of  the  calculations  which  have  already  been  made 
using  the  explicit  scheme. 

(b)  Results  of  Calculations 


To  check  the  accuracy  of  the  program  the  calculations  for  case  No.  1 
(  Vq=  0. 5,  '*'q=  0,  Hq  =  0)  have  been  examined  in  detail.  /  The  pressure  profiles 
at  various  time  intervals  are  depicted  in  Fig.  1,  Since  Tq  =  Mo  =  0  the  Rankine- 
Hugoniot  solution  holds  and,  for  comparison,  these  exact  pressure  values  and 
shock  positions  are  shown  at  t  =  1  and  5  microsec .  It  is  seen  that  the  shock  is 
better  approximated  with  increasing  time,  but  only  2  to  3  microsec.  are  required 
for  sufficient  accuracy. 

The  other  calculated  dependent  variables  may  also  be  compared  with  the 


corresponding  Rankine-Hugoniot  values  for  case  No.  1: 

U 

P 

q 

HUGONIOT 

0.08730 

0.03125 

1.565 

0.2500 

CALCULATED 

0.0867 

0.027 

1.54 

0.25 

The  calculated  quantities  represent  mean  values  about  which  there  are  small 
oscillations  at  the  various  space  mesh  points  behind  the  disturbance.  The  largest 
percentage  discrepancy  is  in  the  specific  internal  energy. 

The  location  where  p  is  hadf  its  majcimum  value  is  considered  to  be  the  shock 
position.  At  each  time  cycle  this  criterion  is  used  to  compute  the  shock  velocity 
in  that  time  increment.  In  Fig.  2  the  calculated  shock  velocities  (relative  to  the 
interface)  for  case  No.  1  are  shown.  As  these  are  found  to  oscillate  widely,  an 
average  shock  velocity  is  computed  according  to  the  formula; 


Shock  Velocity  = 


Position  of  Shock  -  Position  of  Interface 
Time 
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It  is  seen  in  Fig.  2  that  the  shock  velocity  computed  in  this  fashion  converges 
more  rapidly  and  more  smoothly  to  the  Rankine-Hugoniot  solution.  The  ratio  of 
this  value  to  the  true  shock  velocity  for  various  times  is  given  by 

Time  (microsec.)  1  2  3  4  5 

Ratio  1.100  1.034  1.015  1.005  0.999 

In  Figs.  3, 4  the  disturbance  (pulse)  velocities  calculated  in  this  manner  are  depicted 
as  a  function  of  time.  Each  is  seen  to  be  tending  smoothly  towards  an  asymptotic 
value  represented  by  a  dashed  line.  The  value  of  the  asymptote  was  computed 
under  the  rough  assumption  that  the  rate  of  approach  at  corresponding  intervals 
is  the  same  when  |/io|  +  Ko|/0  as  it  was  for  Case  1,  i.e. ,  the  ratios  displayed 
above  for  the  perfect  fluid  calculations  are  employed  in  estimating  all  the  asymp¬ 
totic  values. 


2 

In  Fig.  5  the  value  of  r  is  shown  at  various  time  intervals  for  a  typical 
parameter  combination  (No.  14).  It  is  seen  that  at  each  instant  t2  >To^  only  in  a 
small  region  which  represents  the  current  position  of  the  disturbance.  It  is  only 
in  this  moving  region  of  disturbance  that  the  medium  behaves  as  a  viscous  liquid. 

A  summary  of  all  the  calculations  which  have  been  performed  are  presented 
in  Table  III. 


TABLE  m 

Summary  of  exploratory  calculations  made  with  the  explicit  difference 
scheme,  a  =  1.  5i  A  x  =  0.  1  cm,  At  «  0.  05  |i  sec.  Values  represent 
averages  between  the  interface  and  the  position  of  the  disturbance  after 
5  microseconds.  The  units  are  in  the  gram-centimeter-microsecond 
system. 


No.  1 

^o  ■ 

■  0.  5 

P 

U 

V(.l/p) 

P/Po 

Disturbance 

To 

Ho 

Velocity 

1 

0 

0 

1.  538 

0.0265 

0.0868 

1.466 

0.547 

2 

0 

ao8 

1.  550 

0.0293 

0.  0871 

1.461 

0.549 

9 

ai 

0 

1.  567 

0. 0327 

0. 0876 

1.  452 

0.  560 

10  1 

ai 

ao8 

1,555 

0. 0324 

0.  0877 

1 

1.450 

0.  562 

14  1 

( 

LO 

1  ao8 

1 - i 

1.737 

0.  0681 

0.  0897 

1.418 

0.661 
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IMPLICIT  DIFFERENCE  SCHEME 


In  this  section  an  alternate  computational  scheme  is  outlined  which  is  devised 
to  reduce  the  machine  time  required  for  those  parameter  combinations  with  large 
^^and  Vq  values.  It  is  an  implicit  scheme  based  on  the  Eulerian  formulation, 
whereas  the  original  scheme  is  explicit  and  is  based  on  the  Lagrangian  formula¬ 
tion.  The  new  scheme  is  only  valid  for  impact  between  bodies  of  Identical  material 
with  Therefore  it  does  not  supersede  the  explicit  scheme,  but  serves  as  a 

desirable  complement. 


The  procedure  has  been  adopted  upon  the  suggestion  of  Dr.  Herbert  Keller, 
Institute  of  Mathematical  Sciences,  New  York  University.  In  treating  similar  sys¬ 
tems  of  equations,  Dr.  Keller  has  found  the  implicit  scheme  to  be  unconditionally 
stable.  Thus,  no  restriction  on  At  and  Ax  is  involved,  only  the  desired  accuracy 
need  be  considered  in  choosing  the  increment  sizes. 


(a)  Difference  Equations 


Since  the  two  impacting  bodies  are  identical  the  phenomena  are  symmetric 
about  the  center  of  mass  coordinates  (Fig.  6).  The  calculations  will  be  made  only 
for  bo^  2  where  the  fixed  space  coordinates  are  denoted  by  j  =  0,  1,  2, . . .  There 
dq/dz^O  and,  consequently,  the  Eulerian  formulations  of  the  governing  equations, 
equations  (10)  through  (21)  of  Ref.  2,  reduce  to 


dp 


(3) 

at 

dq 

.  ^  aq 

(4) 

dt 

^^al- 

(5) 

au 

a  t 

^az 

(6) 

p  = 

f  (/»,U) 

A  centered  difference  scheme  is  used  to  provide  more  accuracy.  To  illus¬ 
trate  the  process  the  details  will  be  carried  out  for  equation  (3).  The  scheme  is 
centered  at  point  s  =  j  As,  t  =  (n  +  1/2)  At  as  shown 
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m  ^  —  j 

1 _ 

n  +  1  — ,  1  i 

L  rj 

k  /* 

1 _ 

S  A. 

n  J _ 1 

1  ^ 
1 

7  ' 

: 

«* 

k 

k  /■ 

7 — r- 

1 

1 

1 

\  1  ^ 

h- 

k _ 

L.  . 

U-Az  T 

F  (. 

M  1 

F— 

DISTANCE 


O 

X 


points  used  in  (3),  (5) 
points  used  in  (4) 
denotes  center  point 


in  the  sketch.  Equation  (3)  is  replaced  by 

n  +  1 


i  - 


^■>*1 

J.^1 _ LjLi 

At 


n 


p  —  P 

~i±i _ L±i_  r 

At 


"  I  ' 


-  (/Dq)  "  ■  ^  (/>q)"  -(/oq)*' 


±tl. 


J--  1  +  _ LLL 


= 


2A  z 


2  A  z 


=  0 


Similar  centered  schemes  may  be  written  for  (4)  and  (5).  Upon  simplification 
the  system  becomes 
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(8) 


(9) 


n  +  1  „n  +  l.  n  n 
At  ^^+1  “‘*1-1  ^1  +  1  ^i-1 

"Si  .n  +  1  .  .n  +  1  .  n  .  n 


"vC  ’c 


^j-i  +1  +1  -1 

iFn  +  l.n  +  l  n  ^^n  1  2^orn  +  l  n  +  l.n  n  ll 

"2  h-l*")*!  **’)♦!  + 


The  required  pressure  values  are  computed  from  (6)  aooordlng  to 


(10) 


n  +  l_./n  +  l  ttU  +  In 
Pj-1  **j-i/ 

p?:j  =  ^:}) 

+ 1  ^  + 1’  ^  + 1 ) 
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(b)  Initial  and  Symmetry  Conditions 


The  initial,  n  =  0 ,  values  of  the  dependent  variables  p,  q,  U,  p  are  all  known 
from  the  Ranklne-Hungonlot  equations  (Fig.  7).  To  determine  their  values  at  all 
other  time  intervals  a  method  of  computing  the  values  at  time  n  +  i  from  those 
known  at  time  n  must  be  made  available.  An  iteration  technique  will  be  devised 
for  this  purpose.  ■  The  following  relations,  which  follow  from  the  symmetry  of  the 
problem,  will  be  utilized: 


(11) 


n 

P  .  " 
-] 

n 

"i* 

II 

“  q 

n 

n 

n 

n 

U.  . 

P  ,  " 

P  .  < 

-j 

J  ’ 

-j 

j 

(hence  q  “  =  0 ) 
0 


(c)  Iteration  Procedure 

To  start  the  iteration  let 


(12) 


.s.  n  +  1 

q 

j  (0) 


n 


+  1 

(0) 


n 


n  1  n  1 

be  zero-order  approximations  to  q  |  and  U  "  respectively, 
first-order  approximations  are  calculated  in  the  following  sequence: 


Then  the 
3 


(a-1)  Substitute  ^  ~^j(0)  1*^*°  (^)  *3®°  (1)^  from  the 

resulting  two-term  recurrence  relation  (n  +  1  fixed,  j  varied): 


~n+l^  .s,n+l 
P  *  P 
j-l(l)  j  +  l(l) 


* 

n  ^  n  -A.  r.>,n+l~n+l  ./„\n 

P  ^P  +1^-  ^  q  +v/)q) 

j-1  j  +  1  2  j-l(l)  j-l(o)  j-1 


(13) 


1  4i 

2  AZ 


■.^n+  1 
P 

j  +  l(l) 


.^n  +  1 

q 

j  +  1(0) 


+  (pq) 


n 

j  +  1 


]■ 

« 


3 

This  procedure  is  adopted  upon  the  suggestion  of  Dr.  Herbert  Keller,  Institute  of 
Mathmetical  Sciences,  New  York  University. 
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Several  methods  of  performing  the  calculations  are  discussed  below  in  section  (d). 


(b-1)  To  obtain  a  trial  value  for  ^  use  (10)  with  ^  •  7^ t,}  Md 
„n  +  1  «n  +  1  *  J  J 


(14)  7?  n  +  1  _  n  +  1  «.n  +  1  \ 

(0-1)  To  obtain  a  first-order  approximation  for  use  (9)  wlib  ^ 

p“  ^  replaced  ^  ^  j  ?  j  j  respectively: 


iJ 


n+l 
j-l  0) 


+ 


«n  +  l 
^J  +  i{l) 


I  At  r»vn  +  l  ~n*l  n  _«n+l 

8  Az  [^j-l(o)  +^j  +  l(0)  ^j-KI)^^ 


"  -uP  1 
j+l  Uj-I 


^n  +  i  _^n  +  l 
At  ‘’j  +  l  (6)  ‘’]-l(( 


n  n 


(15)  ■  Az  ~n  +  l  ^  «.n+l  n 

^j-Ml)  ^j  +  IU)'*'  ^j+l  ^j-l 


» 


2x/f  T„ 


1 

"2  [P j- 

+  R 

3  Az  L‘ 


n  +  l  ^n  +  l  ^  n  n  “ 

1(1')  Pj  +  KI')'*’  Pj+l'’'  Pj-I 


n+l 

j+l(0) 


n+l 


n  n  1 

j+r'ij-ij  ■ 
« 


f?n+l 


As  with ^ ,  we  have  here  a  two-term  linear  recurrence  relation  for  *  Its 

solution  is  also  discussed  in  section  (d).  * 

n  +  l 

(d-1)  A  first-order  approximation  for  p  •  is  now  computed  from  (10)  with 

n  +  I  .  ^n-^l  and  ^ 

P  •  -  P  !  in  Ui  -  Wj  (I) 


^  j  (I) 


n+l 


(e-1)  To  calculate  a  first-order  approximation  for  q  •  use  (8)  with 


n  +  l  n+l  n+l  n+l 
— 

for  q  3q/  dz  by  q  [* 

~  n  +  l 

(I)  .  giving 


n  +  l 


^jll)*Pj  "Pj(l)  replace  qj  in  the  difference  expression 
j  (0)  *  latter  substitution  linearizes  the  equation  for 


(17) 


-n"  +l4Lf;rn  +  l  .  „n\ Nn  +  l  ^-n  +  l  n  1 

(1)  8Az  '^j(0)  +1(0)  ”^j-l(0)  ^j  +  l”‘lj-l] 

[pj  +1  (1)  “Pj  -1  (1)  Pj+i  “  Pj  -i] 


lAt  _ 1 

2  Az  <s<n  + 1 


'■'id) 


8  r^yn  +  l  o«a  +  lj.5fa  +  l  .  n  -  n  n 
3  Az  hj+l(l)“2qj  (1) -"^3-1(1) -^qj  +  i-zqj  +qj.i]| 


This  is  a  three -term  linear  recurrence  relation  for  ^  +  1  fixed,  j  varied)* 
a  method  of  solution  is  discussed  below  in  section  (e)/ 


The  second-order  approximations  are  calculated  by  merely  repeating  the 
iteration  process:  (a-2),  (b-2), . . . ,  (e-2).  The  resulting  equations  differ  from  the 
corresponding  equations  of  the  first  iteration  only  in  that  the  subscripts  (1),  (2) 
replace  the  subscripts  (0),  (1)  respectively. 

In  general,  to  proceed  from  the  k-order  approximation  to  the  (k+1)  -  order 
approximation  we  go  through  the  above  iteration  process  (with  subscripts  (0),  (1) 
replaced  by  k  ^d  k+1  respectively).  The  process  is  repeated  until  a  resonable 
convogance  criterion  is  satisfied.  Usually,  only  a  few  cycles  are  required  in 
such  schemes. 


Assume  that  it  has  been  decided  that  K  iterations  are  sufficient.  Then  set 


n  +  l  ~n  +  l  -n  +  l^jn  +  1 
Pj  Pj  (K)’ 


u 

j  (K)’  ^1 


rn  +  l^^n  +  l  „n  +  l^~n  +  l 


1  nn»  Pi 


=  P> 
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and  proceed  to  the  n-t-2  time  step.  A  flow  chart  describing  the  numerical  proce¬ 
dure  is  given  in  Fig.  8. 

(d)  Two -Term  Recurrence  Relations 


It  is  a  consequence  of  the  symmetry  relations  (11)  that  when  J=^0  the  two-term 
recurrence  relation  (a-1)  simplifies  and  yields  the  explicit  formula 


(18) 


~  n  +  1 
'’i(i) 


n 

=  ^1 


1  1  n 


,  .  1  ^t  . 

^■^2  27*1 


mA  -t-  1 
1  (0) 


Now  (a-1)  may  be  used  to  compute  at  all  the  odd  j  space-mesh  points: 

j  (1) 


«<n  +  l 
^  j+2(l) 


r  2  27  *lj+2  (0)J  (1)  ’  2  27  (1)  J 


(19) 


.n  r,  1  At  n  l..n  f,  .  1  At  n1 


Thus,  set  j  =  1  and  compute^  in  terms  of  set  j  =  3  and  compute 


>wn-t-l 


^"■Jj^jintermsof^J-*.  etc. 


#^11  1 

The  p  J  at  even  j  are  determined  by  the  continuity  of  p - for  sufficiently 

large  ^  (the  density  of  the  undisturbed  medium,  which  has  not  yet  been 
reached  by  the  shock  wave).  Let  J  denote  such  a  large  even  integer  at  time 
t  =  (n  + 1)  A  t.  Then  Jij  =  p^  and  the  recursion  formula  is  used  to  calculate 

from  ri^t  to  left.  Thus ,  set  j  =  J  -  2  and  compute  p  ,  ^  /i  \  terms  of  JS^ 

n  +  1  n  +  x""* 

set  j  =  J  -4  and  compute  ?'j_4  in  terms  of  pj  etc. 

A  suitable  value  for  J  may  be  found  by  first  making  the  calculations  for  the  odd 
valued  mesh  points  until  has  decreased  to  the  value  p^.  This  value  is  then 
taken  to  be  J  - 1.  Ihis  left-to-rigdit-to-left  technique  is  depicted  in  the  sketch. 


Calculate  Odd  Stations  X 


j  =9 
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3 

-K- 


Calculate  Even  Stations 


J-4  J-3  J-2  J-1  J 

— M M 1 


■1 
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Alternatively,  the  choice  for  J  may  be  made  from  a  knowledge  of  the  propaga¬ 
tion  velocity  of  the  disturbance.  Then  the  calculations  for  both  the  odd  and  the 
even  stations  may  be  performed  from  right -to-left.  The  advantage  of  this  proce¬ 
dure  is  that  (18)  may  then  be  reserved  to  check  the  calculations. 


The  remarks  made  ^ove  are  also  applicable  to  the  two-term  linear  recur¬ 
rence  relation  for  ,  Here  the  symmetry  relations  together  with  (15)>  with 

J=0,  yield  the  explicit  equation 


(20) 


'jjn  +  1  ^  n  At 
1  (1)  1  Az 


~n  +1  n 
^1  (0) 


yn  +1 
1  (1) 


+  /> 


{-vf 


1  r«{n  +  1  .  ^n  1 .  2  ^o  r»-n  + 
~2  1^1  (IM  Pi  J  3  AT  L 


1 

(0) 


By  replacing  the  word  "density"  by  "apeciflc  internal  energy"  the  remainder  of 
the  dlscuaaion  on  solving  the  recurrence  relation  alao  carries  over. 

(e)  Three -Term  Recurrence  Relation 

Relation  (17)  may  be  rewritten  in  the  form 


(21) 


.  +  1  T.  '^n  +  1 

“■^J  ^  +  1  (1)  ®J  ‘‘j  (1) 


vdiere  Aj,  Bj,  Cj,  Dj  are  known  quantities: 


At 

(Az)2 


.n  +  1 
J  (1) 


+  /» 


n 

j 


(22) 


B, 


1  +  2Aj 


16 


I' 


(22) 


A  method  for  solving  such  a  three-term  linear  reourrenoe  relation  which  is 
particularly  suitable  for  machine  calculations  is  taken  firom  Riohtmyer  (Ref.  6, 

p.  101). 


The  required  inequalities  A.>0>  Bj>0,  C.>0and  Bj>Aj-fCi  are  all  seen  to  be 
satisfied.  The  only  other  requirement  is  that  q  be  speoined  at  the  left  and 
right  hand  boundaries  of  the  space  meah.  J-0  an<r]>J; 


(23) 


-Vo/2 
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CONCLUSIONS 


Preliminary  runs  with  the  program  utilizing  the  explicit  difference  scheme 
have  shown  that  excellent  results  may  be  obtained  provided  the  parameters  Ho 
and  Vq  are  not  too  large.  The  stability  analysis  provides  estimates,  depicted  in 
Table  I,  which  will  serve  as  a  guide  in  the  choice  of  At  for  the  production  runs 
to  be  made  on  the  IBM  7090  at  Eglin  Air  Force  Base.  TTie  stack  of  IBM  cards 
for  this  program  are  being  transmitted  along  with  this  report. 


The  programing  of  the  complementary  implicit  difference  scheme  will  be 
completed.  When  finalized,  it  will  be  sent  to  Eglin  where  production  runs  will 
be  made  for  those  parameter  combinations  with  large  Hq  and  Vq  values.  The 
two  programs  should  permit  the  one-dimensional  computations  for  the  entire 
range  of  impact  parameters  desired. 


Simultaneously,  a  finite  difference  scheme  for  the  visco-plastic  equations 
will  be  sought  for  the  case  of  the  axisymmetrlc  impact  problem.  It  had  been 
anticipated  that  the  particle-in-cell  method  developed  by  Harlow  (Ref.  7)  for 
a  perfect  fluid  could  be  extended  to  our  visco-plastic  equations.  Some  effort 
was  spent  in  this  direction  but  has  currently  been  suspended.  The  method  is 
similar  to  the  one-dimensional  explicit  scheme  in  that  an '’artificial  viscosity'* 
is  introduced  and  it  is  explicit.  Stability  problems  can  therefore  bq  expected 
when  viscosity  is  included.  An  alternative  method  by  Kolsky  (Ref.  8)  may  be 
extended,  but  it  has  given  best  results  when  the  flow  is  nearly  spherical. 


It  may  therefore  be  expedient  to  decompose  the  problem  into  two  parts 
as  suggested  by  the  study  of  the  governing  equations  detailed  in  the  Second 
Quarterly  Report: 

1)  An  initial  stage  in  which  only  the  inertial  terms  of  the  governing  visco¬ 
plastic  model  are  retained.  The  medium  is  then  a  perfect  fluid  and  the 
partlcle-in-cell  method  applies. 

2)  A  secondary  flow  problem  in  which  the  flow  has  smoothed  out  to  a  nearly 
spherical  pattern.  Here  Kolsky's  method  should  yield  good  results. 


The  feasibility  of  such  a  decomposition  will  be  studied  and  a  continued  effort 
will  be  made  to  devise  a  computational  method  valid  during  the  entire  cratering 
process. 
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APPENDIX 


The  explicit  difference  equations  are  given  by  (36)  through  (42)  of  the 
Second  Quarterly  Report.  They  are  as  follows: 
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/n+1  n+1  \ 
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where  we  have  eliminated  X  from  (36)  and  (38)  to  obtain  the  second  equation 
and  the  equation  of  state  is  rewritten  as  p  =  f(V"^,  U)  =  g(V,  U).  Here  -t,  is  a 
parameter  with  dimension  of  length  which  essentially  determines  the  magnitude 
of  the  pseudo-viscosity. 


The  analysis  of  the  stability  of  this  system  follows  the  method  outlined  by 
Richtmyer  (Reference  6).  The  equations  of  firrt  variation  of  (A-1)  will  be 
obtained  in  which  quantities  of  the  second  and  higher  order  are  dropped. 

This  will  give  us  linear  equations  for  the  first  order  variations  4>  S,  (j, 
p  (the  dot  does  not  indicate  time  derivatives)  in  which  the  zero  order  quantities 
appear  as  coefficients.  The  equations  obtained  are 
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.n+1  .n  At  r.n  .n  .  yn  xn  1 

"  Pj^Ax  L^j  +  l/2“^i-l/2'^“j  +  l/2"^j-l/2‘*^j+l/2“®J-l/2j 


At 


t  .n+l  "1 

ST  L^j+i  J 


in+l 
j+1/2  ■  "  3  At 

»+l 


4  4d  1  rxn+1 


i.  r{r«+*  .V"  1 

vLj4l/2  J+I/2J 


<^-">  «^tU=  O  ft"4l‘/2  -^+1/2] 


.n+l 

Vl/2 


(nf  rv‘+*  tv”  1 

jy2  Vst  n’itl/2*''jtl/2j 


2(Po^) 

av 

r^+l 

VAt 

at 

L  j+1/2 

^+1 
j  +  1/2 

au  j+1/2 

The  zero  order  quantities  are  considered  constants  and,  consequently, 
superscripts  and  subscripts  denoting  net  points  are  omitted  from  them. 


The  first  order  quantities  are  asstuned  to  have  the  Fourier  representations 


(A-3) 


— n  ikx 


_n  ikx 


(j?  =Id” 

»  k 

fSl  p”  • 


ikx 

e 


ikx 

e 

ikx 


where  X  =  jAx.  The  von  Neuman  stability  criterion  is,  essentially,  that  the 
coefficients  , .  . . ,  remain  bounded  as  the  calculations  proceed  from 
t^nAt  to  t  =  (n+l)  At.  t  >  (n+2)At,  etc.  To  investigate  this  substitute 
representations  (A-3)  into  (A-2)  and  set  like  harmonics  equal  to  sero  to 
obtain  the  relations 
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(A-4) 

(A-5) 

(A-6) 

(A-7) 

(A-8) 

(A-9) 
where 
(A- 10) 


_n+l  _n  ,  .n+1 

\  =B 


D”"‘=D”-(p+QfS+4|Tj  [b”*'-B"] 


,n+l 


(Po^)' 

ZV^ 


/ivf  ■ 

1  2(p^i) 

1  r  ° 

av 

r“1 

\at  -)  l\  . 

J  ^  V  At 

at 

L^k 

k  av  k  a  u  k 


8  = 


-  2  8in(kAx/2) 


p  Ax 
o 


Relations  (A-5),  (A-6)  may  be  combined  to  yield 
-n  4  -  n 

This  may  be  substituted  into  (A-4)  to  obtain 
(A-ll)  Allf  ‘  =  (l  -i  ^  [e”  ,  F”  ] 

If  (A-ll)  is  substituted  into  (A*i)the  result  is 
(A-12)  +  1  P  At  ^1  —  3  At)  Aj^+  3  (At)  J 

Now  (A-12)  may  be  used  to  eliminate  from  (A-7)  and  (A-8)  with  the  results 

Iv 
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»n+l 


(A- 14) 


2V^ 

2(p^t) 


V  A  t 


(Irf  (*  -I  T 


Finally,  substitution  of  (A- 12),  (A- 14)  into  (A-9)  yields 

4, 


(A- 15) 


Upon  introduction  of  the  notations 


g  ■  2  sin  (kAx/2) 

0  Ax 
o 


(A<*  16) 


A  »  p+Q+S+i/j 


e=-i-iy 

®  2v  at 


li  - 

A^ 

av 

^au 

(Po^) 

2 

av 

V 

”  at 

the  relations  (A- 11)  through  (A-15)  may  be  written  in  the  matrix  form 


(A- 17) 
where 


v“*‘  .  gy” 


- 1 

_ 1 

,n+l 

\ 

n+1 

®k 

Y®'*’'  m 

n+1 

k 

< 

-n+l 

*^k 

^-k 

n+1 

'k 

...  — 
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and  the  "amplification  matrix",  which  depends  on  and  T^,  is  given  by 


(A-18) 


Y 

iPV^t 

-ipVA^t 

ipYCAt(-|^-e) 

ipVAAt 


■ipAt 


0  -ipAt 
0  ^(At)^ 

1  -^A(At)^ 

-2eC  0  ^c(At)^(^ -e)  0^C(At)^(^  -e) 


e^At)^ 

-P^A(At)^ 


Sg  5g 

hv  au 


e^6(At)^ 


^6(At>' 


Now,  expansion  of  the  determinant  |G  -  Xl|  and  setting  the  result  equal 
to  zero  shows  that  the  eigenvalues,  X,  of  G  satisfy  the  equation 

(A-19)  X^(X-l)|(X-l)^-(X-l)^At[6At  +  C(2-€At)-|  +a^(At)^2eC-6)|  •  0 

2  2 

1£  t~  constant  and  At/(Ax)  s  0(1}  as  At,  Ax^^O,  then  P  At  ■*  0(1)  and  the 
secular  equation  reduces  to 

(A-20)  X^X-l)^|x-l-P^At[2C-|  *Y]-0(At)J 

Thus,  the  von  Neuman  requirement  for  stability 
(A-21) 

is  satisfied  provided 


X^  I  s  1+0  (At) 
m  I 


£  1 


where  we  have  used  the  fact  that  dV/at  <  0.  The  inequality  ^11  hold  provided 
At 


(A-22) 


(Ax)^  |^^+4(p^t)^|||| 


2idV| 


It  may  be  noted  from  (A-22)  that  in  the  limit  as  the  viscosity  tends  to 
zero  the  stability  criterion  reduces  to 


(A-23) 


At 


(Ax) 


which  is  the  value  given  by  Richtmyer  (Ref.  6,  p.  220)  for  this  case. 
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In  the  case  ^  ^  0  it  is  also  possible  to  derive  a  stability  criterion  from  (A-19) 
and  (A-21)  under  the  condition  that  a  Ax,  where  a  is  some  constant,  instead 
of  holding  t  constant.  Then  the  restriction  is  found  to  be  relaxed  to 


(A-24) 


^  V^l|v  ■  ^  0  lu' 


Criterion  (A-23)  holds  in  the  region  of  the  shock,  and  criterion  (A-24)  is  valid 
in  regions  away  from  the  shock.  However,  if  true  viscosity  is  present  (|^  ^  0) 
then  we  must  always  have  At/(Ax)^  «  0(1)  for  stability;  if  At/ Ax  *  0(1)  one  of 
the  eigenvalues  goes  to  infinity.  This  results  since  in  this  case  the  ^  term 
in  (A- 19)  does  not  have  a  factor  of  order  (Ax)^  but  .of  order  unity. 
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FIG.  I  CALCULATED  PRESSURE  PROFILE  Ax  «O.I.  At»0.05y 

COMFW^O  TO  THE  RANKINE-HUGONIOT  SOLUTION 
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WITH  ASYMPTOTIC  VALUES 


CO 
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o 

H 


CO 
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LU 
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CO  CD 

tr 


to  CO 
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UJ 
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CL 
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u. 


(b) 

F1G.6  ILLUSTRATION  OF  IMPACT  SITUATION  SHOWING 
(a)CENTER  OF  MASS  COORDINATES,(b)EULERIAN 

SPACE-MESH  POINTS 
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BODY  1 


BODY  2 


Station 


-  00  to  -1 

0 

0 

1  to  ® 

Vo/2 

0 

0 

-  Vo/2 

Po 

Pod-Po/d)"^ 

Po 

0 

Vo^/8 

Vo^/S 

0 

0 

(M 

0 

> 

v„2d/4 

D 

ere  d  is  determined  by  the  relation  (Ref,  Z,  eq,  63) 


d  «  £( . -  ) 

1-Po/d  ^ 


Fig.  7  Tabular  display  of  initial  and  boundary  data  for  implicit 

difference  scheme. 


FIG.  8  FLOW  CHART  DEPICTING  THE  SEQUENCE  OF 
CALCULATIONS  IN  THE  IMPLICIT  DIFFERENCE  SCHEME 


APGC-TN-61-43 


INITIAL  DISTRIBUTION 


1  NASA 

3  OAR  (RROSA) 

3  ASD(ASOR) 

1  AEC(Tech  Info  Serv) 

2  US  Nav  Rsch  Lab 

2  US  Nav  Ord  Lab 

1  ASD(ASRMDS-21) 

2  BSD(BSRVE) 

2  SSD(SSRTA-l) 

2  Ballistic  Rsch  Lab 

10  ASTIA(TIPCR) 

APGC 
3  ASQW 

2  PGAPI 

3  PGEH 

20  PGWR 


Page  34  of  34  pages 


«  0  U  U  «  ± 

«  if  2  S  1 S  -a 

X  K  <  U  Q.  £  0. 

-  ->ri  ^>>5; 


^  to  («  O  7i 

t  i-  “•  “5  £  .!i  S  3 

Is^  uiiii: 

-  -  B  ^  >  p 


I  u 

=  -  Se 

i  S  s.  “i  -  & 

§5«-  “ 

.  o  qS  ^  ^  T, 

ir  {14  >.507.2 

o  0  73  2 

rt  Z  «  *7  )S  5 
o  0  !3  •  :£.  « 
•2  G  1  3  g  u 

•  5  h  ^  < 

"  S  '■2  "  ■ 
EiS^  5l 
•a  >  -1 3  t 

m  Z  U  3  S 

4,  Q  <« 

S  jS  g  £ 
a,i5  “  § 
wy  y  “-U 
*  0^  C  S 
>*  W  ^  S* 

B  •«  <0  O 

1=33  s 

“  -s  o  ..  1 

u  S  u  «  ? 

•r  >  (0  ^  t. 

sis  §“• 

« -2 
■3  O'O  ?i- 

"  Z2“iz 

0.  Q  H  N  H 


5  o  a  ..  .0  5 

B  ^  'B  ^ 

c  ^  2  S  t)  C  >0 

5  o  A  2  u  •*  fl 

Uu*"^(j*c  0* 

«  U  41  ft.  C  c 

6  41  „  5  3  5  O 

1  -a  -s :  S 1 1  s 
1^ il * » ®  • 
“!§shit 

5222  0S.3, 

>«o-ft.S*>*-K- 

®  •  e  «>  'O  B  #®  5 

§  s ! :  s  a‘“:2 

8^3.5 

p  B  II  4|  *<  c  «• 

®5 

a  «  •  2  «  ^ 

si-IS|;= 

S  Sfii  3-3 
asPs^jSS 

"3  'S  j  8  0  «  £ 

•S  a  a  t  J  S  S  3 

-S“8.?S-a5 

2  S  -  2  0  S  •  • 

AB>a,ur4lh 

tisSiin 

siJliSlI 


'  >  V  >  .ti 

E •-  it 

3  «  g  I* 

^  X  rt  a 

'S  •  ft  'S 

o  j3  ■'<  S  'B 
S--E  “  >- 
£  u  -a  *3 

j  "'  2  3  B 

B  5  » I  2 

W  «  ^  'B  P 

§  e  9  B  3 
J;  .44  ^  y 

•3!5o«- 

or  u  A  ^  .e 

4t  «  o  h  u 

±  w  2 

.2  S  *«  *  > 

5  52^6 

•  •  *0  .44  l« 

2  §  s  'a  >• 

ia  i 

l,i  •  *■ 

•  ffl  a  a 

ll-SS? 

^  £  -u  g 

•H  41  ‘O  -• 

«  1 1  "  S 

S  a^ 

246*1 

is-l?  s 

m3  '  e  > 

. "  3  •  o  -a 
2  S  c  ■£  • 

S  .M  4>  *0  '*■ 
^5  h  e  «4 
.?  _  Jl  3  *4 

SIIH 


$  o 

»  -  St 
?  ■  S 

2  W  “a!  _  £• 

»  z  U  Q  "  “ 

^O  0!  J,  .  TJ 

Sotsl? 

tj  |i  •I.'  *A  ‘2 

2  5  t  •  2  • 

o  O  «  JS  ^  jS 

•B  r  5  3  <»  o 

sgll< 

"  2  *  • 


Sis*^  ;l 

o  >•  j  3  t 

CO  Z  U)  -H  u 

„  ”  Q  -  14 

s  js  a  £ 

tB  s  o  ®‘0 

>.  W  H 

ii^-i 

|z5:s!: 

S  o  11  2 

u  5  O  V  ? 

£p|£ 
s  !<^.2 

£“:2  Sx? 

I)  M  to  a  S  T 
c  .2  2  A  ?  j, 
n  ■  W  •  7  'B 

^  “so;.' 
b  M  .  2  z 

Ou  Q  H  N  H 


£  2  3  5  V 

0  i  S  5  B 

s  3:  1-ji 

o  “  fl  “  2  ! 

.2  SS  2  »• 

c  «  «  o  B  B 

II  h  a " 


®  5  B  4»  TJ  B 
B  5  4)  •  C  2 
0  9  t)  .  (4  5 

•fl  -.  <B  g  _  b 

42  :'S :  s 

1 1  3  -S  s  b 

b  ^  4*  4  5  b4 

£  U  I  B  •  TJ 

^  5  g  •  2  *> 
s :  •£  g  2  a 

2  8  e  •s  -° ! 
S  !  2  -  £ 

■3  f  .  8  0  • 

•>  «  b  ?J  5 

■i  S  i  i 

a  2  4»  0  ■ 

*  'S  e< ^  Jl , 

"3  8  ^  S  ?  *4 

pltll 


o  >*  8 

^2  i  5  I 
2  9  J  B  a 
5  (<  u  a  4>  ^ 

"  5  b  „  X  (4 

®  .  H  frjq  J: 

'  •  «  B  u  rs 

t>  b  ,  •3  2  5 

a  S  q  s'®  I 

'•6525 

b  .  u  ”*  -B  •*“ 
«  b  b  ♦;  .  *0 

M  3  '3  S  O  4* 

S  8  S'S  a2 

.3  •  V  4>  Ob 

b  -  b  b 

0  s  ■  4)  I,  r 

S  S2'S 

.  0  , ,  ,  • 

3.5 -Mil 

2  S  •  ?|  6 

■o  B  u  M  •  2 

2  ®  I  5  *  S 

5  g  t  S  X 
9  O  *  b  U 

838  II 2 
«£ i -  LS 
£12^6* 

4  B  I.  ..  B  . 

.  B  *.  l.'S  = 

5  *  .s  3  B  s 

»  -‘  •  S  c  2 

••  » 'S  5  B  B 

:  girti  8 

5  P  j  ,•  5  o 

5  .A  8  5  9  § 


■§“|  <s 

M  Cm  $  *4  ^  * 

s  ^  o  U  O  -  • 

>  .bCO^fUtlo 

h  >*Q.U.  b  «v.5 

S  41m.**  '9* 
flueb^Boab 
>4  -S  .M  B  b  C  <« 
rc  £  ^  o  ft.  £  ft. 


"  a  ri  >  >  5: 


f  .  b  (0  (t  W  •*  y 

S  5-  ft.  Ib  b  •  W  .S 
8.qi.<eOa,B 

>s  2  2  0  B  e  4 

X  S  <  u  0.  2  £ 

--a  B>;>:c 


z  o 

»  -  St 

?  ■  S  8. 

5  2  Q  2  * 

•  O  A  ^  ^  .Q 

b  ib  2*15  7  2 

2012^3 

8  z  t  ^3  4 

:2b8< 

S  i-  2  4  . 

Syis 

S  J  O  I  £ 

ar5=t '  8 

w  3  u  “'O 

a2!S  -I 

c  2  ^  A  ®' 

£  Z  Qa  ^  ^ 

-a  fi  " 

u  S  U  «  o 
r  3  to  A  b 

1  —  c?  0  (1 
u  >■  ^  s 

511*2 

!“:2  sr5 

c  .2  z  2  «  * 

2  5  W  .  r-  3 


212  S  t 

'M®Ll 

I-*! 

risSia 
1-38 :  *1 

7| 1 8  E : 


I  b  -  C  • 
o  0  M  M  .44  M 

b  b  b  0  « 

■s  ^  " 

®  2  0  ®  *0  .B, 
a  0  «  •  B  3 
0  P  4>  <4  S 

•d-4-B  g-.  b 

tmii 

Elsas'® 

A  E  C  "  >0 
i*  4»  «  3  • 

s  :«§.£■§ 

•  fr'g-.s 

«  0  Jb  *3  3  ^ 

a  J{  H  g  ^  2 

^  4:  3  S  - 

4»  -  b  tl  5 

■a|8.S^i 
-  1  «  S  .?  £ 

'G  8  .S'  0  0  tJ 

ll  8  1'32 

b  .  b  4*  'g  b 


>«.?  >  4. 

^  E  '3 

e  A  9  ,B  I 

S  <•  u  <a  4 

o5.£.s5 

1  «  g  I 
b.^03- 

o  h  a.j 

•  •  b  £  4 
b  b  .  —  3 

.2  5  S  3I 

g  •  S'  ♦ 

^  £  .S  i 

b  b,  •  •• 

•  b  b  - 

's*is; 

"abb 

.  .*  b  b  . 

0  «  •  b  I 

3  si 
_:  ®  •  2 
3  2-S.|  ! 

-■‘Mi 

2  2  •  ?  i 

8  §2  ?.: 
:  B I  a  .4 

9  O  b  b  t 
g  'fl  .0  2  j 

•  V  «  A  1 
:fl  •-  B  O  8 


b  w 

2  1®.“ 

•38.3 

4  .  ■■'  U 
b  b  ■  tfl 

b  g  '2  !2 

B  i  ?.'3 


«2q"T'?  STbb.j^bS^’s 

ls“3^ 


S  w 
3 

?.«?£» 

5  •  Q  3  £ 
•  o  «.  >  r-  4 
JTib  >•«  7  2 
2  O  B  2  -3 

2  z  r  •  £  • 

So  S 

sMU 

m  z  u  a  s 

!  JO  I  £ 

w  a  u  “"U 

.  0!  5  B.  — 

>^u  H  "  — 

5S2  S 

li^is 
I  ^  ?  3 
“ISti 

a  i  <  £ 

B  .2  £  2  «  ^ 
5  J  w  .  »7  3 

U  M  «■  n  I 


5  S£  .. 

r-  g  • 

f  1 2  S  i 

t  s7  at E 
S  ■*  I  . .“  £ 

1.88:81 
ii  ii  I  o 

« B  a  s  E  B 

o  •  M  *«  ■**  4fi 

5  g  o  B  B 

S  b  h  b  0  g 

a  5.  g  >°  b 

®  •  a  g  “g  .5 

B  0  g  •  B  3 

0  P  b  .  (4  3 

■rt  _  A  i!  b 

S  b  M  .2  •  4» 

•8  'S  i  s  3  8. 

g  g  A  a  b  B 

g  g  g  o  -e  '3 

g  ■s.  • «  *2^  1 

B  1  2  4  0  0 

B  ®  js  a  ^ 

a  g  3  j  a 

«  -f  .  ■3  g  S 
<1  _  t  s  • 
■6 1  a  g  1 3 

l3  ®  B  0  ^  4 
.B  C  ^  B  .B  B 

g  g  A*  o  n  •* 

2  £  b  'X  y  g 

M  g  0  -o  g 


•  GTbg-gi; 

£i“=5  ilWil 


>®3  S  3  P 
■o^i  B  -a  0 

g  A  9  B  B  Q 
M  (4  g  a  g  ^ 

0=  .£  a  «  * 
*<“248 
^  0  3  •**  S 

®  1 

'  ■  V  E  ^4  '3 

ih'iil 

8 si:i: 

b  w  g 

g  b  b  *;  .V 

'  f  S  i  8  A  •* 

•3  •  g  g  e  b 

b  -  *•  b  A  « 

,®  S  2  2  b  g 

3'|l.|1| 

2S  •  f|  6 
•o  e  u  M  «  g 
g  o  .5  fl  .2  £ 

Isfi-. 

o  *5  ■«  g  3  b 

£  :2| 5  • 

a  ^  4  3  S 

2  .£  u  V  0 

*  B  •  S  8  3 
’”  >  •  a  £  1 
:*£•’£  8 

llSIiil 


irt  is  given  for  the  machine  program  which  is  currently  being  written.  I  I  chart  is  given  for  the  machine  program  which  is  currently  being  written. 


4  r  O 

u  Q  ?  s 
O  .  .5  ■<  *  . 

M  Li4  5  M  9^  ■ 

s  0  O  u  „  • 

t.-£ss,rs-s 

s  S  ^  o  £  I  iS 


i  u 

fl  '  J!  «  5.  2 
.2  u  2  cd  ^  2* 

k..*  S  0«  “  2?  ti 

*  2  «  Q  h 

^o  »!>"•>. 

2  2  r  "•  2  ! 

■2  H  1 3  g  o 

•  a  M  ■?  < 

V  ^  -H  5 

sa^-o 

o  >■  ►^  3  Z 

tf5  2  W  •«<  u 

sjS’S  S 
g-S^-s  § 
w  y  u 

lll-s 

c  5  r  ^ 

g  ^  P<  ^  *1 

u  .S  o  "  • 


U  2  U  V  O 

||g|£ 

^“■a  S-" 

s  i  s  2  g  ’ 

E  2  °  •»  z 

0.  3  H  3  H 


J  s  ..  ii- 

■  S  *1^3 

eg  *0  'S 

S  S  «  d  >2 
3  2  X  •«  3 
^  a  -i  p  o“ 
5  r  ~ 

a  5  '§  a.  °  • 
a  “  o  a  •  « 

O  V  4)  V  h 

3  °  >  .C  9  « 
U  •  ■*3  *2  > 

lls^ss 

a  g  E !  i  s 

?  2  -  5  S  ^ 
u  i*  o  2  -^  ■ 

X  «  >  h  ^ 

B  «  *0  .S  ,®  s 
u  ■  c  2  bi  ^ 
u  .  4  iS  «M 

8  fc  .0 

.•2:5s* 

3  • ! "  a  .a 

ajill! 

5^  fti-S 

.  3  ?  S  5  I 

c  “  a  •  •  s 

•  Its i  i 


3  a  o  g  S  », 
I  &  •  a  S 
•  **'$««« 

mill 


>  «l  »  -H- 

£  -s  0  > 

3  ,g  c  a  * 
o  *a  «  ^  ^ 

.2  " 

3  •  d  I 

g  g  3  a 
•  •-  2  2  0 
c  2  £ 

«  ^  2 
S  £  2  5  H 

4j  •»•  .O  ••«  U 

^  H  4  *0  « 

•3So^* 
w  o  A  r3  .d 
g  g  o  ** 

^  u 

•  2  g  g  > 

mis 

jl^lr 

g  g  g 

2  ?.:^2 

I  0  X  ^  ^ 
^  h  g  g 

:||i2s 

4  S  4  >.  - 

S.?“s  8 

a  3 :  i  i 

'.:9  81  “ 

lllil 


•  ^  S  u  S  -  • 

*  ►££  2  8  s  a 

^ah^sodah 

z  2  '<3  6  £  £  £ 
aSi>? 


w  >*  Z  n 

o  -«  .  a;  - 

?  •  I i 

■Is  S' 17  • 

.  O  Hi  ^  ^  .o 

4  <4  ‘•s  7  a 

sSs-s’s 

s  8  S aS-a 

•  saT-si? 
s  t;  2  «  . 
8i2^  ;S 

■S  >  J  2  “ 

ui  z  u  S  S 

5  JO  5  s 

aj  s|g 
w  a  o  “'o 

.  otf  44  *■  •* 

&SS"f 

I?  3  3  “ 

B  2  (1«  ^  *4 

O  ^  u 

-s  o  4  .* 

•p  g  H  ? 
si?  §“• 

a  S  j(Sl^ 

«  SS  -2 

3  o  *  ^  ST 

S  i  I  'I  S 

a  !  «  -73 
22«2z 

Q.  Q  H  N  H 


s  *0  2  £  g  d  <2 
a  §  A  S  ^  *  5 

4  a  “  !■;  E  i<>* 

5  ■*  a  •  .s  2  ~  ^ 

g  g  -o  0  5  a  8  8 

1-8  11  "3  g  8 
S'!  3  S  f  8  •  ? 


mniii 

2  •  y*!*  —  5  5  ^ 

•4  »■  a.  S  >  *•  4  •’ 

^■Qai«d02 

aggS2*3U,<* 

o  9  V  .  M  3  •« 

•a  _  .a  4  _  4  .0 

8  15 

£®4»j|^a-0- 

5*^6"  —  "Oflx 

g*iiBsgg**Od' 

UOVO'C’JVO. 

g  *3.  •  S  S|  •  ^ 

J;q22  0*03.® 
JJgCS^^wXJ' 

a  8  u  j  5  a  * 
■3'Sj3o"’®i 

s  1  5  >•  S I  S  !  ■ 

cga.ta3a4 

A2v2>k."«g 

Sc  2--‘-'3® 
g  g  ^  o  •  -  2  .  • 

3..lt'3.*8sg. 


E  s  *  I 

3  ,£  c  a  * 

g  A  g  ^  M 
.£-£*•? 
!3  •  ?  1 

fr-fi  Jl 
g  u  z  a 
.  -5  2  2  ? 

®  2  » I  £ 

I  a  i  a  3 

3 :  o  •  -2 

O'  g  A  Z  X 
g  g  A  k  g 

s  2  h  *  » 
5  -'S’g 
•  2  2  "  £ 
■Mil? 
’S'  3  a  •  ^ 

•?18S 

|?‘ei 

I  St  g 

8  >  I  a  • 

a  JS  E  •  5 
t  a  4  >,*■ 


g  g  ^  O  z  *4  5  . 

Il|!|lp 

g  I  *•  *  >  g  V  «> 

osSa"".** 


k  ■»  «  >4^ 

a1  |4 

•sisla 

M  -  g  o  k 

Sliil 


4.  >-2  3 

o  ■“  .  a  4 
? , s  ?i 0 

Is  8.“=^®' 

§  ^  Q  -  “ 

*  O  ^  >  t**  “O 

k  tk  ►•‘2  7  £ 

SSis^iS 

2  2  £  T  ^  m 
o  0  5  js  ^  jS 
•5  C  3  X  •  g 

2ga|-< 

sa^ 

g  >  J  S  t 

to  2  U  g 

,,  ■-'  Q  «  It 

£  J  o  g  5 

s.5=^i§ 

.  etf  «  'g' 

=  5«"f 

liSdl 

4 1 8  s  •? 

VJt- 

4  §  4  i 

8  .3  i  2  «  ' 

n  S  “  -73 

£  0  h  M  H 


d  g 

-  >«.?  >  «  > 

•  3  i  -i .  I  s 

£S3«««®bA 

^  t!  ^  4j  X  g  d 

B  0  *  •-  .A  »  _  .S 
gd.^Wg.'gg 

2v4^0Z*Sx 

s  ®  .  ^  frx  j:  >. 

I*  •  g  S  .£  •-  Z 
t,  g  k  .•'<236 

S-3  8  S  S'^l  2 

SS4*<4^K 

“’•B.giss 

.  k  w  g  •'*  X  —  g 
•g^k**  .‘0« 
g  M  I;  -A  •  «*  4>  " 

c2>B'gAtSx 
<«Zagg  Akg 
k^.Akk^g-g 

a  ,0  a  •  »  4  4  ■» 

g  z  ®  •  5  2  *  2 
U  '2  S  -S  8  g  2  M 
h  2  .£  4  -3  2  I  o 
Svarste'^Xk 

-•^'g  B'3>‘g  «- 
oexg^sPg 

g-AA  ggCg 

I  8  8^  S-S  E 

^  g  A  •«  *  Z  iJ 
Vaug>'S2£ 

aSk  £  **,4^ 

•  gae*>gsa 

2kgAxs  2 

5  4  k  k  •*  J  >4^ 

•  a  g  S  ►‘•g  3  S 

g  ■  5  a  5  £ 

*'‘g^a3^s.A 

C  g  J  •  a  S  3  “ 

5  «  ►  ’5.5  £  I  « 

Slisliil 


■£  lI  *4  ^  * 

**^0044  4 

>  .kioagt'y 

s  t“'  S  s  •  3  ■= 

P>>Sk<B0S.k 

z  2  '2  ij  £  I  £ 


i  o 

^  -  2- 
?  •  S  fi  2 
Is  ^°‘s» 

*  2  g  Q  2 

•  0  Ag^r-o 

4  a-  a-S  7  a 
a  o  4  S  —2 

_  •’*’!£  a 

|§5aia 

l%hli 

zhj< 

|a";| 

w  Z  U  =  u 

8  j  8  g  2 

3.5='-:  i 

«  a  u  “-o 

4  A  ^ 

S-SS"? 

I  2  fc  S'  *4 

^  -s  o :  8 

g  S  U  g  o 
-r  g  w  X  k 

t  •£  5  2  tt. 

8  •  j  S  o 

aM  C  ^  W  (A 

U  e  <  £ 

alSg  i<~ 

8  .3  I  '2  fT 
a  !  «  -73 
“  850^' 

£  S  H  3  H 


S12  S  S 

S  e  A  k  V 
8  8  g  Its 

8  ■*  I  *  .?  S 

|.st:g| 

illll: 

•  g  2  ►  g  g 

i>  a'  •A  w 

•9  k  r  k  o  4 

k*  0.  5.  g  k 
®  •  B  g  *0  .£ 

I*lsi2 

pj-sr 

-a  I  4  2  8 

g  4  u  o  -e  '3 

g  i  X  -g  0  o 

■3'?j3o- 

•S  s  I,ga3 

®  'S  «l  s  t  * 

•A  B  4  k  -2  ► 

u  g  ^  0  !  «4 

llllsi 

8T  4  ••S'S 

ism& 


i.  >•  g  s 

>®.A  ►  2  ► 

^  Z  E  *3  Ok 

c  9  2  c  c  ^ 

"aaa?l-| 

^<o“la|l 

»  .  H  a...  *  >. 

a  g  B  «  -a  Z 
g  k  .  2  fl  fl 

2  S  6  a  I  •  S 

g  i  I  "  '  t 

►  a  g  B  3  B  3 
t,  ^  5  -  X  -  u 
^  u  ^  .‘da 

?g3S®il* 

3  ►  9*  V  A  :S  X 
3  •  a  g  Akg 

^  .2  k  k  4®«  I 

®  Z  •  V  k  g  9 

3  sa^  g 

.  a  .  4  q  •  g 
|S's|1?  a- 

u  m  ^  ^  JS  L. 


lUim 

^2  ^  8 1  -s  ^ 

S  4  "  4  2  ^ 
8  85|  |S| 

a  g  ^  ^  o  3  g 

vHiH 

srsisit 

:  sir'!  s  g  « 

I  Ja  Sis  II 


